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Non-equilibrium fluctuations of various stochastic variables, such as work and entropy production, 
have been widely discussed recently in the context of large deviations, cumulants and fluctuation 
relations. Typically, one looks at the distribution of these observables, at large fixed time. To 
characterize the precise stochastic nature of the process, we here address the distribution in the 
time domain. In particular, we focus on the first passage time distribution (FPTD) of entropy 
production, in several realistic models. We find that the fluctuation relation symmetry plays a crucial 
role in getting the typical asymptotic behavior. Similarities and differences to the simple random 
walk picture are discussed. For a driven particle in the ring geometry, the mean residence time is 
connected to the particle current and the steady state distribution, and it leads to a fluctuation 
relation-like symmetry in terms of the FPTD. 

PACS numbers: 05.40.-a,05.40. Jc,05.70.Ln 


Introduction. — The past two decades have witnessed 
significant development in nonequilibrium thermody¬ 
namics M. The fluctuation relations are remarkable 
discoveries which have quantitatively refined the concept 
of the second law as applied to small systems. One 
of the central issues in nonequilibrium statistical physics 
has been in characterizing the universal nature of fluctua¬ 
tions of thermodynamic variables, such as heat and work 
that quantify entropy generated in non-equilibrium pro¬ 
cesses. Usually, one measures the accumulated entropic 
variable, say X , over a fixed time interval r, and its fluc¬ 
tuations are then characterized through a distribution 
P{X). Defining, for example, X as the stochastic total 
entropy, one can then prove the detailed and integral type 
of fluctuation relation for any fixed time interval r, in var¬ 
ious Markov processes ;9j. For large observation times, 
one finds the large deviation form P(X) ~ e Th ( x / T ) 0, 
where h(x) is the large deviation function. The corre¬ 
sponding cumulant generating function (CGF) is defined 
by /x(£) = log {e^ x )/r where (...) is an average over the 
steady state, and this generates the nth order of cumu¬ 
lant I n as 


In = 9"KC)/rk=0. 


(1) 


For physical quantities related to entropy production, it 
is well known that the CGF shows the fluctuation re¬ 


lation symmetry [ill [12]. This symmetry is not only 
mathematically beautiful but also physically important 
since it reproduces linear response results and al so give s 
nontrivial relationships on nonlinear responses |13H15|. 
Large deviations and the CGF have been crucial towards 
constructing universal thermod yna mic structure of the 
nonequilibrium steady state 16-18|. 

The large deviation function gives us the probability 
of observing rare events in some fixed observation time 
window. An interesting and natural question to ask is as 
to how long would one have to wait to see a rare event of 


a specified size?. This is just the question of the first pas¬ 
sage problem for the stochastic variable X. Although the 
physics of fluctuation at fixed time has been intensively 
studied and a lot of discoveries have been made, sur¬ 
prisingly, only a little is known on the stochastic nature 
of its time evolution itself. One expects that the time- 
evolution of stochastic thermodynamic variables should 
behave like a biased random walk in some configuration 
space, but the details of the temporal aspects have not 
been investigated. 


The main aim of this Letter is to investigate this as¬ 
pect, which is clearly necessary for a deeper understand¬ 
ing of stochastic thermodynamics. In particular we con¬ 
sider the problem of the first passage time distribution 
(FPTD) of the desired stochastic variable, which is an 
experimentally measurable quantity. The FPTD here is 
the distribution of waiting time at which a stochastic 
variable first reaches some target value. We consider the 
typical properties of the FPTD of entropy-related vari¬ 
ables within the broad and well-established paradigms of 
stochastic thermodynamics. Three examples of nonequi¬ 
librium processes are considered: (a) an over-damped 
driven particle in a ring geometry, (b) classical charge 
transfer via a quantum-dot, and (c) heat transfer across 
a coupled oscillator system [see Fig. m for a schematic 
description]. Note that due to recent development in 
time-resolved measurement techniques, there are a num¬ 
ber of relevant experiments for these setups that look at 
nonequilibrium fluctuations lfil- 241. 


Using these models, we address the following questions. 
Is there a typical functional form for the FPTD and espe¬ 
cially its tail ? How does it depend on the sign of entropy 
produced ? What are the differences as compared to the 
time-evolution of a simple biased random walk ? Con¬ 
cerning this last question, consider the case where X is 
the position of a biased diffusing particle on the open 
line. Then, defining J- rw (t, X)dt as the probability that 
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FIG. 1: (color online) Schematic picture of setups, (a): The 
over-damped driven particle in the ring geometry. The right 
picture shows the potential landscape in an infinite line pic¬ 
ture. (b): Classical charge transfer via a quantum-dot. (c): 
Heat transfer via a coupled oscillator system. 


the particle hits X for the first time between times t and 
t + dt , one easily finds jUJ 


Frw(t,X) = 


\X\e 


-(X-ht) 2 /2I 2 t 


i°g* 


\Z2nl2t 3 


( 2 ) 


realization of the stochastic process, the particle makes 
the transition (zq 0) — > ( a,N ) in time t. Then the work 
done is w = /(NL + a — v) while the heat dissipated into 
the bath is q = w — U a + U v . Thus, at sufficiently large 
times, entropy production rate is proportional to the av¬ 
erage rate of the winding number. Due to the positive 
force /, the particle on average moves in the positive di¬ 
rection and the average winding number rate is positive. 
However there is a finite probability to observe the parti¬ 
cle moving in the opposite direction. The ratio of proba¬ 
bilities between positive and negative winding number at 
any finite time is quantitatively given by the fluctuation 
relation. We address the FPTD for the winding number. 
Let T a ^(N,t) be the transition probability from (y, 0) 
to (a, TV) and let F a ^(N,t) be the probability that first 
passage between (z/, 0) to (a, N) occurs between time t 
to t + dt. We note the relation for N ^ 0 

Ta :V (N, t) = f du T ajOl (0,t — u)F a:l/ (TV, u). (4) 

Jo 

Taking the Laplace transformation T a ^(N,s) = 
/ 0 °° dt e~ st T a ^(N,t) and similarly for the FPTD we get 


We will use this as a reference form, and aim to figure 
out the similarities and dissimilarities between this sim¬ 
ple random walk picture, and the real stochastic time- 
evolution of entropic variables. Intriguingly with use of 
the fluctuation relation symmetry, one can derive the 
asymptotic form of the FPTD for these models, and 
dissimilarities to F rw can be argued. In addition, we 
derive the exact expression for mean residence time for 
the driven particle in the ring geometry which leads to 
fluctuation-relation-like equality in terms of the FPTD. 

Driven particle in the ring geometry. — We consider a 
colloidal particle driven by a constant force and confined 
to move on a ring, as depicted in Fig. (la). The dynamics 
is well-described by the over-damped Langevin equations 
with temperature /3 _1 . The Boltzmann constant is set to 
unity and let us also set / > 0. To proceed, we discretize 
the space into L sites on the ring separated by a small 
spacing a. Let P v {t) be the probability to find the parti¬ 
cle on the zzth site at time t. Its evolution is given by 

dP (t) 

—^r 1 = W v v -\P v -i(t) + W„, v+1 P u+1 (t) - W VtV P v {t ), 
at 

( 3 ) 

where W Vy0l is the transition rate matrix element 
which satisfies the local detailed balance condition 
Wv+i^/Wv^+i = e -Pi u v+-L- u v)+P a f an d u v i s the po¬ 
tential energy at the zAh site. It is useful to introduce 
the winding number, N, which, for any given particle tra¬ 
jectory, is obtained by counting the number of times the 
particle makes the transition from site L to the first site 
(reverse transitions from the first site to the site L are 
counted with a negative sign). The particle’s state can be 
labeled by the duplet (z/, N). Suppose that in any given 


F a , v {N,t) 


1 

27T* 


/»c+ioo 


ds e s 


Tg,u{N, s) 
T a , a (0,s) ' 


( 5 ) 


We consider the asymptotic behavior of the FPTD at 
sufficiently large waiting time. To this end, one can write 
the time-evolution equation of the joint probability of 
the variables zq N. Solving this through Fourier-Laplace 
transformation, one can get the formal expression for the 
transition probability matrix j26| 


T(JV '*> = 2S 


dz 


A(z 


z N+1 det [s — W z ] 


( 6 ) 


where the matrix W z is given by the matrix W replacing 
(1,L) and (L, 1) elements by zWi^l and z _1 Wl ,i respec¬ 
tively, and A(z, s) is the co-factor matrix for the matrix 
of denominator. There are two singular values z±(s) from 
the denominator, which are connected to each other by 
the fluctuation relation symmetry 261 

( 7 ) 


2+(s)z-(s) = e J fLa , 


where fdfLa is the entropy produced in the reservoir for 
a single winding around the ring. Setting a = zq using 
the above symmetry and Eq. one can express the 
distribution in terms of only one singular point [26j] 

F„ „(iV, t) = / ds e^ s+b log2+(s)I , (8) 

2m J 

where b = N/t. The steadty state FPTD of winding 
is then given by F(N,t) = Fv,v{N, t) P„ S , with the 
steady state distribution p^ s . A further careful exami¬ 
nation reveals that the singular value z+ (s) is connected 
to the CGF, p.(£), for the winding number 

s — n(£) = 0 , where £ = log 2 + (s) . (9) 
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Based on these relations, a saddle point analysis leads to 
the following exact asymptotic expression of the FPTD 
oc J r asym (t) for sufficiently large waiting time [2(3 


-Fasym (t) = exp[-Ff - (3/2) log t] , 

°° (-T, \n+2 

r = E^T2 )T^)I C =o (10) 

If I 3 I* (31$ - hh)lf 

2 1 2 61 f 24/f 


where the function g„(£) is connected to the CGF; 

„ (p\— d ] n 

— [\ d £ 2 ) d{] V d I s > ' 

Some important observations on Eq. m are now in 
order. The asymptotic temporal decay form depends 
neither on the sign nor on the amplitude of the wind¬ 
ing number, although the actual probability of negative 
and positive winding numbers differ by exponential fac¬ 
tor (in entropy produced). Thus, even extremely rare 
events follow the same asymptotic form. In the linear re¬ 
sponse regime with small first cumulant, the asymptotic 
behavior is well-explained by the simple random walk 
picture T rw . In the far-from-equilibrium regime, how¬ 
ever, critical deviation from this picture reveals itself in 
the higher order terms with nontrivial expressions. This 
deviation will be significant in small systems where the 
degree of nonequilibrium is easily increased. We note 
that the asymptotic form is given by the general form, 
in terms of cumulants, irrespective of detailed potential 
forms. This indicates that it might be applicable to wider 
classes of physical situations. As we see below, it turns 
out that the expression is valid for many other situations 
when the cumulants are calculated for appropriate phys¬ 
ical quantities. 

Two other examples. — We now show that the asymp¬ 
totic form m also appears for open nonequilibrium 
systems such as (b) classical charge transport via a 
quantum-dot and (c) heat transfer via coupled oscilla¬ 
tors [See the Figs. (lb,lc) for schematic pictures]. 

Case (b): Let pl and p R be respectively the chemical 
potential for the left and right leads and consider spin-less 
electrons transmitted via a quantum-dot with an onsite 
energy e. We measure transmitted electron at the right 
contact to the reservoir, and let the accumulated electron 
transfer till time t be n. Charge transfer produces Joule 
heating and is directly connected to entropy production 
rate as (S) = f)(p L - p R )(h). 

We now consider the FPTD of the accumulation of 
electron number, an experimentally measurable quantity. 
Let 1 and 2 respectively denote the unoccupied and occu¬ 
pied states of the quantum-dot. Then the time-evolution 
of the two states is described by the same type of dy¬ 
namics as in Eq. ©. The transition probability Wqj is 
composed of two contributions from the left and right 
reservoirs W = W L + W R , where W{ 2 = 7 [1 — / r (e)] 
and Wf i = "/fr(e) where f r is the Fermi distribution of 


the rth lead (r = L,R). Hence these elements satisfy the 
detailed balance LFf 2/^21 = The modified 

transition probability matrix W z in (|6| is given by the 
matrix W on replacing Wf\ and W R X in the off-diagonal 
matrix elements by W R 2 z and W R iZ _1 respectively. The 
singular points in the denominator are £±(s) which are 
again connected by the fluctuation relation symmetry 

^_(s)z+(s) = (11) 

In the present example it is easy to see that the first 
passage from the initial state (i, n = 0) to any fixed 
desired value of n, also fixes the final configuration j. 
Using the renewal equation, we can obtain FPTD from 
(*, 0) —> (j, n), using the same argument as for the driven 
particle in the ring geometry, and find that the FPTD 
is proportional to Eq. m where now the cumulants are 
for charge transfer and known exactly (see j26j). For 
the case of many sites with strong onsite-Coulomb inter¬ 
action, one may employ the symmetric simple exclusion 
process [ 2 ^]. We can demonstrate that the same expres¬ 
sion is obtained analytically for this system of coupled 
quantum-dots. 

Case(c): We consider the example of the coupled os¬ 
cillator system, exchanging heat with two heat reservoirs 
at temperatures T Rl T Rl whose dynamics is described by 
the overdamped Langevin equation 

7 x 1 = —kx + r) L (t) , 7±2 = kx + rj R (t) , (12) 


where x \ t 2 are the positions of the first and second par¬ 
ticles which are coupled via spring constant k , and x = 
X\ — X 2 - The noise terms r] r satisfy the fluctuation dis¬ 
sipation relations (Vr(t)v r '(t')) = 25r-,r'7/3 j r 1 J(t — t’). In 
this case we consider the heat transfer into the right bath 
in time t and this is given by Q = f 0 dt'kx(t')[kx(t') + 
r lR(t')\/'y and are interested in the FPTD for transition 
from an initial state (x, Q = 0) to a state with Q amount 
of heat transferred. As in the previous examples, we can 
think of our system executing biased diffusion in (x, Q) 
space. However in this case fixing Q an the initial x does 
not fix the final position and so an extension of the for¬ 
mulation is required. A heuristic derivation is given in 
[i(3 | but the final result for the tail of the FPTD for Q 
turns out to be the same as given by Eq. Co| where now 
the cumulants for heat transfer are known exactly from 

Uli 


Numerical demonstration of the asymptotic formula for 
several cases. — We numerically verify the asymptotic be¬ 
havior ca for the three examples discussed and shown 
in Fig. 0- In the case (a), we consider the dynamics in 
continuous space given by jx = —dU(x)/dx + f + pit), 
where we impose the periodic boundary condition with 
the periodic length 1 and we employ the potential U(x) = 
sin( 47 ra;)/( 47 r). The variable r/(t) is the Langevin noise 
satisfying {f](t)v(t')) = ‘2jP~ 1 5(t — t'). In case (b), we nu¬ 
merically update the state using a Monte-Carlo approach 
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FIG. 2: (color online) The FPTDs for the three models are 
shown with unit normalization. In the inset in (a), unnormal¬ 
ized data are also shown, which shows that the FPTD with 
the negative entropy production is very small. The black solid 
lines are fitting by the theory m, while the dotted line is 
the asymptotic curve of the random walk ©• Parameters 
sets are (a): j3 = 5.0, (b): {pl,Pr) = (6.0, 1.0) and (c): 
( Pl,/3r ) = (5.0,1.0). All other parameters are set to one. 


with the specified transition rates. For (c) the system 
evolves through the Langevin dynamics in Eq. (1121) . In 
all cases we sample the initial state from the steady state 
and then measure the FPTD for specified values of wind¬ 
ing number [in case (a)], the charge transfer into right 
reservoir [in case (b)] and the heat transfer [in case (c)]. 

For events with a negative entropy production there is 
a finite probability of the event not occurring at all in 
a given realization. Hence for negative values, we plot 
the distribution, conditioned on the probability that it 
occurs. The results are shown in the Fig. (O. In short 
time scale, non-universal behavior is observed. However, 
all three cases clearly show that the asymptotic behavior 
is well-described by the theory m, irrespective of the 
fixed values including even very rare events. At finite 
times the logarithmic correction is important. The devi¬ 
ation from the simple random walk picture (which gives 
T = lf/(2I 2 )) is also clear. 

Basic equation and integral fluctuation relation in 
terms of first passage. — Let us consider the entropy pro¬ 
duced in the thermal reservoir for case (a). Let T a>ll (S, t ) 
be the transition probability from the site zz to a in time 
t during which the entropy of the heat reservoir increases 
by the amount S. The entropy is determined by the po¬ 
tential at the sites a and v and the work done by exter¬ 
nal force, i.e., S = ft [U v —U a + fa(a — v + LN)}. Note 
that fixing zz, S does not necessarily fix a, N. Let us de¬ 


fine lFa,v{S,t) to be the FPTD only for S [not for (a, £>)], 
while reaching d from v. For fixed 5, the site a depends 
on zz. It is uniquely determined if |<S| is sufficiently large, 
while for small |<S|, there can be at most two choices of 
a, respectively on the two sides of zz. As example, see the 
right figure in Fig.([TJi) of a case where two a (denoted by 
a ±) can be reached for a fixed negative <S, starting from 
the site zz. Then we note the following basic equation in 
the Laplace representation 

Ta, v (S,s) = y^T a ,a'(0,s) J 7 q/, l/ (5, s) ■ (13) 


This type of equations provides in general, a basis for 
considering the FPTD for entropic variables. 

We now establish several relations. The first is an exact 
relation for the mean residence time at a given lattice 
point for given entropy production, given by 


/ 0 °° dtT ajU (-S,t) = e s p^ s /J, 
/ 0 °° dt T^ a (S, t) = p^ s / J , 


where p^ s is the steady state distribution at the site a 
and J is the steady state particle current. In the first 
relation, it is assumed that the process (zz, 0) —> (a, —S) 
is opposite to the direction of current, while in the sec¬ 
ond relation (a,0) —> (zz, S) is in the same direction as 
current. These are connected by the detailed fluctuation 
relation [27|. Note that the sign of S is not specified. 
The proof for this is presented in [26j]. Eas. (fl3l) and (fldl) 
are crucial for deriving other relations as we now show. 

We now employ the usual definition of total entropy 
= ln{p^ s /p‘^ s ) +N • Then for fixed negative entropy 
S < 0, using relations (fT3l) and (THl) leads to the equality 
J2a e ~ Sa ’ u -^a,u(S,s = 0) = 1. Multiplying both sides of 
this equation by p^ s and summing over zz immediately 
leads to the integral type of fluctuation relation 


«e 


»s = 1, 


(15) 


where the average ((■■■))s implies taking all possible first 
passage paths producing the negative entropy S, and that 
start from the steady state. Numerical demonstration is 
presented in [26]. 

Summary. — In general, it is difficult to characterize 
general temporal aspects in stochastic time-evolution of 
thermally fluctuating objects. As a first step in this di¬ 
rection, we consider the first passage time distribution 
of entropy-production in several models that are relevant 
to recent experimental setups. We find the asymptotic 
behavior of Eq. m, which seems to be the typical func¬ 
tional form, valid in many situations. For the paradig¬ 
matic example of a particle driven round a periodic po¬ 
tential, we find further properties, given by (fT4l) and (fldl) . 
that characterize the FPTD. It is proposed that Ea. (fl3l) 
is in general the basic equation needed while considering 
the FPTD for entropic variables. 
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Derivation of Eq. (6) in the main text .— Let us define 
the winding number A of a typical trajectory of the dif¬ 
fusing particle by counting the number of times it makes 
the transition from the site L to the first site, L —> 1. 
The winding number decreases for the reverse transition 
1 —> L. Let P„(N,t) be the joint probability that the 
particle is at the i/ th site with winding number N at time 


t. We are interested in finding the probability vector 

P(N,t ) = ,P L (N,t)} .(16) 

Then it is easy to see that this joint probability satisfies 
the following equation 


dP(N,t) 

Ft 


W-P(N - 1, t) + W+P{N + 1 ,t) + W 0 P(N , t ). 


(17) 


Here W is a L x L matrix whose only non-vanishing 
element is L = Wi,l, W+ is a L x L matrix 

whose only non-vanishing element is [W+k, = W L , i, 
and W 0 = W - W- - W+. 

We define the generating function 

OO r. 

P(z,s ) = / dt z N e~ st P(N,t). (18) 

N=-oo ** 

Then one readily find that this satisfies the equality 

P{z,s) = s Po (19) 

where W z = W-Z + W+z- 1 + W 0 and po is the initial 
condition po = P(N = 0, t = 0). From this, one gets 


Overall structure and fluctuation relation symmetry 

In the expression ||221) . crucial roles are played by the 
singular points in the denominator. We first note that 
the determinant has the functional form 

det [s — W z ] = - nfc=1 Wk+hk [z - *+(*)] [z - z-(s )], 

(23) 

where z + (s) and Z-(s) are the singular points. By look¬ 
ing at the tridigonal matrix W z , one finds the relation 
between these singular points 

z+(s)z-(s) = = e -hf L * . (24) 

1 1 Wk+l ’ k 


P(N,t) = 


P(N,s ) = 


/»c+ioo 


2iri 

1 


dse ts P(N , s), 

200 

dz 1 

■Po • 


2iri T z N+1 s-W. 


Hence the transition matrix T(N , s) is given by 


T(N,s) = 


1 

27 xi 
1 


dz 


1 


Z N + 1 g _ ^ 

dz A(z,s ) 


27ri J z N+1 det [s — W z ] 


( 20 ) 

( 21 ) 


( 22 ) 


where the matrix A{z, s ) is the cofactor matrix of s—W z . 
Derivation of the FPTD, Eq. (10), in the main text .— 


Hence we can take these points as |z+(s)| > 1 and 
(s)| < 1. In the limit of s —> 0, z+(s) —> 1. The 
relation (l24l) corresponds to the fluctuation relation in 
terms of winding number which stands for the entropy 
PfLa generated for every increase of the winding num¬ 
ber. 

We consider the matrix element of the transition ma¬ 
trix T(N,s). We first note that the cofactor matrix has 
dependence on either z or z° or z -1 . Hence a matrix 
element is given by the following type of integration 

27 Tij >dz [z-z + (s)\[z-z-{s)\ 1 
where n = IV — 1, IV, IV + 1 and C is a constant depen- 
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dent on the matrix element. We here note 

I dz z~ n (s) _ ( ’ n ^° 

/ 2tu [2-2+(s)] [ z - 2_( s )] I -="W „ <n 

(26) 

This implies that depending the sign of the winding num¬ 
ber N (| AT | > 1), the expression of numerator takes either 
z+ or However, using the symmetry (1241) we can in 
unified way express those only with the singular point 

z +: 


eigenvalue of the matrix W z . Thus if A k(z) is the k th 
eigenvalue of W z and Afc = o is, say, the largest eigenvalue 
then we have 

Mf) = A 0 (z) , with £ = log 2 . (34) 

The eigenvalues Xk are given by the roots of the de- 
terminantal eqation det [A — W~\ = 0. This is identi¬ 
cal to the equation for finding the roots z + , 2 _, namely 
det [s — W z \ = 0 if we replacing A by s. We also note 
that Ao(z) —► 0 as 2 —► 1. Using this continuity in terms 
of the variable s around s = 0, we see that the singular 
value z+ is related to the CGF via the relation 


z _|n| ( s') 

T a ,v(N , s) = C a ,„(a) + , (27) 

2 _(s) — 2 +(s) 

where the prefactor C a _ v (s) accounts for the amplitude 
of transition. For instance, between positive an negative 
winding number there is exponentially large difference in 
the amplitude of the prefactor. 


Saddle point analysis 


S - M (£) = 0,where £ = log z + (s*(b )). (35) 

Thus the function h(b) is completely specified by the fol¬ 
lowing equations 


h(b) = b£(s*) - s * , 

s* = MO * 

7 I 


ds 


- 1 = 0 . 


(36) 

(37) 

(38) 


Using the relation [Eq (5) in main text] between the 
FPTD and the transition probability, we then get 

CvW) = A./*e» 0.(28) 

Using Eq. (EH), we discuss the asymptotic behavior of the 
FPTD. At large t. the function F„ tl/ (N,t) is given by 

F v , v [N,t) °c Jdse-M'V , (29) 

g{s,b) = -s + b\ogz+(s), (30) 

where b = \n\/t. We make the saddle point analysis, 
where the saddle point s* satisfies 

_ 1 b dlogz+(s) 
ds 



From the last two equations it is easy to see that b = 
dfi/d £ and hence it is clear that h(b) is the LDF corre¬ 
sponding to the CGF MO- 

We now express the value T = h(b = 0) in terms of 
physical quantities. To this end, we expand the func¬ 
tion in a Taylor series around its maximum, b m , which 
satisfies 


dh(b)/db 


b=bm 


log z+{s*(b m )) = 0 . (39) 


This implies z+(s(b m )) = 1, hence s = 0 and therefore 
h{bm) = 0. Also £ = 0 and b = d/i/d £ implies that 
b m = Iii the first cumulant of the winding number. Thus 
we get 


T = h(b = 0) = J2 (-0) fe M fc) (/i)/fc!. (40) 

k—2 


where :A_(s) = dz+(s)/ds. As shown in the next subsec¬ 
tion g”{s*) oc 1/b 2 . From this, one gets 

F v>v (N,t) -> e -tUfc=0)-( 3 /2)io g t, (32) 

h(b) = g(s*(b), b) = —s*(b) + blogz + (s*(b)), (33) 

where we introduced the function h{b) to emphasize that 
the saddle point s* is a function of b. We noted in Ea. (l3^1) 
that asymptotic behavior implies b —> 0. 

The function h(b) is in fact precisely the large deviation 
function (LDF). To see this, we first note that the cumu¬ 
lant generating function (CGF) is given by the largest 


The final task is to express the derivatives in 

terms of cumulants of the winding number. To derive 
the expression of the second derivative M 2 )(/i), we start 
with the expression h^ 2 \b m ) = 6 _1 ds* /d&| b=bm . Using 
the relation b = d/i/dt; one gets ds*/d&| b=bm = h/h- 
Hence 

h {2 \h) = y. (41) 

Higher order terms are systematically derived in a similar 
manner and we get 
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hW(b m ) = l/I 2 (42) 

h {3) (b m ) = -h/I 3 (43) 

h^(b m ) = [3 / 3 2 - I 2 I A ] /if (44) 

h {5 \b m ) = — [l5lf — 10 / 2 / 3/4 + if I 5 ] /if (45) 

h (6 \b m ) = [l05/ 3 4 - 105I 2 I 3 2 I 4 + 10lfl 4 2 + 15/|/ 3 / 5 - l!h\ III (46) 


r 


In conclusion we get the asymptotic behavior of the 
FPTD F asym (t ) 

F asym (t) = e —r *—( 3 / 2 ) log *, (47) 

r = £ (n + 2 )! > ( 48 ) 

where <?„(£) = (^^) -1 ^ (^r^) -1 - In the main 

text, we wrote the expression with up to l/ 4 / 

Logarithmic correction term: g"{s*) oc 1/& 2 

We consider the equation (l23l) for determining z+(s). 
We consider the structure of the quadratic equation 


z 2 +c{s)z + d = 0 , 


(49) 


where d = fl/Li W k , k+ i/W k+ i,k = e pfLa . Then the 
solution is 


z+(s) = — c(s) + \Jc 2 (s) — 4 d / 2 . 


(50) 


Hence the function g(s,b) is given by 
g(s,b) = s + &log| 


-c(s) + \/cP(sj~4d /2|.(51) 


The first derivative is then given by 
dg(s,b) 


ds 


= 1 + 6 - 


dc(s ) 

-| I c(s) 

ds 

yj c 2 (s)— 4d 


2 z+{s) 


(52) 


Now we consider the case of b -C 1. For the above to 
be zero, the term ^ycF(s)~^4d must be extremely small. 
Hence, we make the rough estimate 


y/c 2 {s*) -4doc b. (53) 

The second derivative is then estimated to be 


d 2 g(s,b) 


ds 2 


oc 6 [c 2 (s*) — 4d] 32 ocl/ 6 2 . (54) 


Derivation of Eq. (If) in the main text .— We first 
note that the steady state distribution and current can be 
exactly solved for the driven particle in the ring geometry. 


vl s = 


1 


1 "Wot,oi— l^^a — l,a:—2 


Wa- 2,a-l W a -2,a-lW a -3,ot-2 

L L 

j = n w a , a+1 \/z, 


Ot—l 


ot—1 


(55) 

(56) 


where we used the notation Wij = Wi+L,j+L and Z terms of steady state and currents. 


is the normalization factor. There are two approaches 


In the Laplace representation for the time-domain, 


towards getting the expression of mean residence time in P(TV, s), the Ea. (fl7l) is reduced to 


sP{N,s) = W-P(N-1 ,s) + W + P(N+1,s) + W 0 P(N,s) + 6 Ni oPo. 


(57) 


For N 0, let us of the form 

P(N,s ) = z~ N V, 


( 58 ) 


where V is a constant vector. Plugging this into (1571) for 
TV ^ 0 gives the following equation for determining z and 
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V 


[s~W z ]V = 0. (59) 


A careful look at these equations reveals that there are 
two sets of solutions to these equations. To see this, we 
write the above equation in the following form: 


f 3+[W L ,1 + W 2 ,i] -Z+ 

^ -zi u 


where 



0 . (60) 


Z+ = (Wi.2,0 ,--- ,zWi,l), (61) 

Z_ = (W 2 ,i,0,--- ,«- 1 W l , 1 ) t , (62) 

U = sI L -i — W 2L , (63) 

V' = (V 2 ,V 3r -- ,V l ) t . (64) 


Here the matrix W 2L denotes (L— 1) x (L— 1) sub-matrix 
of W excluding the first row and column, while Il-i is 
unit matrix of dimension (L — 1). We have set the first 
element Vi to one. Then we get the following equations 
for V and z 


8 + W l , 1 + W2 ,i-W 1 ,2V 2 -zW 1 , l V l = 0, (65) 

V' = U~ l Z_. (66) 


The second equation leads to the relation V a +i = 
[U~?W 2t i + U~l_ 1 z~ 1 W L .i\ , for a = 1 , ■ ■ ■ , L- 1 . Since 
U does not depend on z we see that V a are linear 
functions of z _1 . Hence putting back V 2 ,Vl into the 
first equation above, we get a quadratic equation for 
z. For the two solutions we get two corresponding ex¬ 
plicit forms for the vectors V. We denote the two solu¬ 
tions by (z + (s), V + (s)} and (z_(s), V“(s)}. From the 
equation for z, we see the fluctuation relation symmetry 
z+( S )z_( S ) = nLi W k , k+1 /W k+1 , k = e~Pf La (< 1). 

Let us now look for a solution corresponding to the 
initial condition that the point starts from v with N = 0. 
A possible solution of Ea. (l57l) is 


f A+z+ n V+ for N > 0 
P(N, s) = < A_zZ N V~ for iV < 0 (67) 

F° for N = 0. 


The unknown constants A + , A_ and V° can fixed by re¬ 
quiring that our above solution satisfies Ea. (l57l) at the 
sites corresponding to v and its two nearest neighbors. 
Clearly then the vector V° must have the following struc¬ 
ture 


F° = (A-Vr,A-V 2 ~,..- ,A-V-_ 1 ,A 0 ,A+V+ +1 ,.~ ,A + V+) . (68) 


There are three constants (A_, A + , Ao) to be determined 
and these will follow by writing the three special equa¬ 
tions at the site v and its neighbors. Let us assume, for 


the moment, that none of these three sites is a boundary 
site on the cell (i.e., v — 1 > 1, i/ + 1 < L). Then we get 
the following equations by looking at the block of N = 0 


+ = W v - 1>v Ao, (69) 

[a + Wv-i, v + W U+ 1 , V \ Ao - W^-iV-^A- - W„, v+1 V+ +1 A + = 1, (70) 

{[s + W v ,v+x + W v + 2,i/+l] v+ +1 ~ Wu+ i, v+ 2 V+ +2 } A+ = W V+1 , V A 0 , (71) 


Using the equation satisfied by which is given from 
the block of N ^ 0, we find that the first and third 
equations yield 


Plugging these into the middle equation, one gets A 0 : 


An = 


s T Wi/~ T W v +i,u 

Vi; 


v,v+V 


K + +i 


W + 


= w. 


-1 

iy,iy—l 


rK 


v+l 


Kiu - 1 


L Vu 


V+ 


A- — Aq/V v , A + — A 0 /V+ . 


(72) 


(73) 
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Special case s = 0. For this case the equation W Z V = 0 
has two solutions. One deary is for z+ = 1 and this is 
the steady state solution so we can choose 

V+ = p ss . (74) 

The other solution for — n*=i Wk, , k+ i/W k+1>k is 
given by 

/ W 2 ,1 W 2 , 1 W 3 , 2 W 2A ---Wl, l . 1 \ t 

" 1 , W lt2 , W h2 W 2 y"' , W lt2 ---W L _ hL ) A ' 

as can be easily verified. From (T73l) we then get 


Tlo 




r Wu+i^ 

■ W v , v+ 1 


Kid - 1 




(76) 


Using the fact that J = [W„+ i, v V+ — W VtV+ iKh_ 1 ] we 
then get 


4> = V+/J. (77) 

From Ea. d72l) we get 

A_ = V+/JV~, A+ = l/J, (78) 


From this we finally obtain, for the case s = 0, the fol¬ 
lowing transition matrices for any states ct, v where v is 
the one “down the hill” (i.e. current is in the direction 
a —> v). 


T a , a {N = 0 , s = 0 ) 
T a ,v{N = 0 , s = 0 ) 
T Via (N = 0,s = 0) 



We finally explain how to obtain Eq. (14), in the main 
text, using Eq. (USD. We note that the entropy produced 
in the thermal reservoir for the process a —> v is given by 
S = /3 [U a — U v + fa(v — a + LN )]. This implies that 
the process (a, N = 0) —> (a,N = 0) is equivalent to 


the process (a, S = 0) — > ( a,S = 0). Hence Ea. d79l) is 
equivalent to 

T a , a {S = 0,3 = 0) = ^- . (82) 

Now consider the process (a, S = 0) —► (iq<S) whose 
direction is the same as the average current. In this case 
we note F V}Cl (S, s = 0) = 1 which means that the process 
will occur with probaility one. This gives 


T„ t a(S,s = 0) = F v>a (S, s = 0)T„ >v (S = 0, s = 0) 


vl s 

J 


(83) 


In the backward process (zq <S = ())—► (a, —S) which is 
opposite to the direction of average current, the detailed 
fluctuation relation immediately gives 


T a , v (-S,s = 0 ) 


,-sPi 


ss 


J 


(84) 


These give Eq. (14) in the main text. 

The FPTD of charge transfer via a quantum-dot .— 
The dynamics of the charge transfer in classical transport 
is described by 


dP(t)/dt = WP(t ), (85) 

W = W r , (86) 

r=L,R 

where both W and W r are 2x2 matrices. P\ (t) and 
P 2 (t) are respectively stand for the probability for un¬ 
occupied and occupied state inside the quantum-dot. 
Standard setup takes the transition matrix element as 
Wf 2 = 7 [1 — /r(e)] and y = 7 f r (e) where f r (e) is 
the Fermi-distribution for the rth reservoir and 7 is a 
hopping rate. Hence it satisfies the detailed balance 
Wf 2/^27 = e /3 ^- /ir ^. Without loss of generality, we 
can impose hl > /r r . 

Let Pi(Q,t) be the joint probability for the i (= 1 , 2 ) 
state (=unoccupied or occupied state) and transmitted 
charge Q measured at the right reservoirs till time t. The 
dynamics is given by 


d t Pi{Q,f) = W 1 , 1 P 1 (Q,t) + W^ 2 P 2 (Q 1 t) + W^ 2 P 2 (Q-l,t), (87) 

d t P 2 (Q,t) = W 2 , 2 P 2 (Q,t) + W 21 Pi(Q 1 t) + bF 2 fl 1 Pi(Q + !,*)■ ( 88 ) 


We define the generating function 

00 pOO 

P(z,s) = J2 / dtz^e- st P{Q,t), (89) 

Q=-oo“'° 


From the dynamics for the joint probabilities, this is for¬ 
mally given by the expression 

p ( z ,s) = s _ 1 w Po , (90) 
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where W z is given by 


W z = 


W n 

W 2 l i + W 2 r \ z~ 


W{; 2 


w< 


W R 2 z 


22 


The transition matrix is hence given by 
T(Q,s) = 


(91) 


1 q 

[ dZ 1 

(92) 

2m J 

Z Q+ 1 s-W z ’ 

1 „ 

f dz A(z,s) 

(93) 

2m ] 

z Q+i d e t _ W~\ 


Now one can see the same structure to the case of ring 
geometry. From this, we find the singularities ~±(s) by 
solving the equation det [s — W z ] = 0, which are con¬ 
nected by the fluctuation relation symmetry 


z-(s)z+(s) 


W£ 2 W* i 

W r 2 W 2 l i 


— g -P(»L-Vr) 


(94) 


In the same way as in the driven particle in the ring 
geometry, any matrix elements of the FPTD have the 
following dependence 


^ pc—ioo 

F a ,„(Q,t) = — / dse st T a>v (Q,s)/T a!<x (0,s) 

J c+ioo 


OC 


27 xi 


ds e ~ t9{s ' b) 


g (s,b) = s b log Z-\. (s ). 


(95) 

(96) 


where b = \Q\/t and C is time-independent matrix. The 
argument from this point follows the calculations from 
Eas. (l29l) and (151)1) in the driven particle in the ring ge¬ 
ometry. Hence, we can reach the same expression as in 
Eas. CTl) and C51) . 

Heuristic derivation of the FPTD asymptotic .— Con¬ 
sider a general process with discrete configuration space 
X and let us look at the joint distribution of X and 
some quantity Q (like heat or charge). This distribution, 
P(X, Q, t), will satisfy the equation of motion 


dP(X,Q,t) 

dt 


C Q P(X,Q,t) , 


(97) 


while the generating function Z(X,£) 
f dQe^P(X, Q) will satisfy the equation 


=CzZ(X,t,t) . (98) 

The large time solution of this equation with the initial 
condition X = Y,<3 = 0atf = 0is given by 


T x , Y (fy t) ~ e^V(X, 0 X (Y, 0 , (99) 


where p, <j>, x are respectively the largest eigenvalue of C j, 
and the corresponding right and left eigenvectors. Taking 
a time-Laplace transformation we get 


#x,fl x(Y,0 

[s-m 


( 100 ) 


Taking an inverse Laplace transformation in the variable 
Q, we get 

Tx.y(Q,s) ~ e-« + W'«l </»(X,e + (s)) x(Y,e+(s)) (101) 

where we assume, based on empirical observations, that 
only the singularity, which satisfies the following re¬ 
lation, contributes: 


p(t)-s = 0. (102) 

We now consider first passage only of the variable Q with¬ 
out caring for the configuration coordinates X. Let us de¬ 
fine Jx,y (Q,t) as the probability that the system starts 
from Y at time t = 0 with Q = 0, first reaches Q in 
the time interval (t, t + dt) and is at X during that time 
interval. Then we have, for Q fy 0, 

P-x,y(Q,s) = ^Tx,x'(Q = 0,s)Txc Y (Q,s) • (103) 

X' 


Now using Eq. (film and assuming that the wave- 
functions do not contribute to the asymptotic behaviour 
we get 


^x. Y (Q,s)~e-« + ^l Q l . (104) 

Finally, transforming back to the time domain, and doing 
a saddle-point analysis, we get 

J r x lY (0,*)~e-®W t , 

where g(b) = £ + (s*)b — s* , b=Q/t , (105) 

and s* is determined by 

s* = p(0, ^7 5 ~ 1 = 0 - ( 106 ) 

Thus we recover the required Eqs. (1361371381) . 

We here note that in all models (a)-(c) the equation 
(I102D yields the quadratic equation of the following type 

z 2 (s) + c(s)z(s) + d = 0, (107) 

which gives two solutions z±(s) and z(s) is connected 
to £(s) by the relation £(s) = logz(s), and the constant 
term d comes from the fluctuation relation symmetry (as 
in Ea. (H51) for the case of driven particle in the ring ge¬ 
ometry). Then from the same argument as in Sec.II-C, 
the same logarithmic correction is obtained. Hence, we 
get asymptotic behavior fill) and (1481) . 

Integral fluctuation relation in terms of first passage 
and numerical verification .— We consider the entropy 
produced in the thermal reservoir S for the driven par¬ 
ticle in the ring geometry. Let T & ,v{S,t) be the FPTD 
only for S [not for (a, 5)], while reaching a from v. De¬ 
pending on (y 1 S), there are two possible situations; in the 
first case a is unique, and in the second case there are 


?x,Y(fys) 
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two values a, as depicted in the right figure in Fig. (la) Hence, in any cases, we have the identity 
in the main text. In both these cases we note the relation 


T a ^(S,s) = ^T S;a /(0,s) T a ', v {S, s). (108) 

ot' 


Then we consider ]T) g J r a , j!/ (5, s = 0) for negative 

S. We first consider the case where a is unique. Note 
that in this case (i/,S = 0) —> (a,S) is opposite to the 
direction of the average current. Then 


^Te- S “ 7^(5,0) = 1 . 


This immediately leads to 


(115) 


e- s ftr &iV (S,s = 0) 


e 


-sPa S T a ^(S,s = 0) 
P s „ s T a , a (0,s = 0) 


n ss p s n ss 
sPa e Py 


rtS S 


P& S 


1, (109) 


where we used Eqs. (|5^]) - ([M|) to get the final result. 

We next consider the case where a has two choices. 
We call a± these two points, located respectively on the 
positive and negative side of v [See the right figure in 
Fig. (la) in the main text]. Using Eq. (11081) . we have the 
following relation 


( Ta+,v («S> 0) \ / 7 a+i [ / (iS, 0) \ 

\T s _ tV (S, 0) ) \T a _^(S,0) J ’ 

where the matrix M is given by 


( 110 ) 


( T a+ , a+ { 0,0) T a+ , a _(0,0)\ 
V 7 a _, a+ (0,0) T a _, a _(0,0) J 


(111) 


We note the following expressions in terms of the steady 
state distribution and steady state current 


( T s+ , V (S, 0 ) \ 

\ 7 a _ !iy ( l S, 0) J 

M 


1 

J 

1 

J 


p ss 

pfe s 


S S ss 

ra.* ra-L. 

SS ss 

"a- 


( 112 ) 

(113) 


Using these expressions, straight forward calculation 
leads to 


«e- St0t )) S = 1, (5 < 0). (116) 


We finally present the numerical demonstration of 
(11161) as well as the relation Eq. ([82]) . The Langevin 
equation was numerically solved with the same parame¬ 
ters set as in Fig. (2a) in the main text. 



0 0.25 0.5 0.75 1 

x 

(color online) Numerical demonstration of Ea. (l82l) in 
the main graph and fluctuation-relation-like symmetry 
in the inset. The Langevin equation was numerically 
solved with the same parameters set as in Fig. (2a) in 
the main text. 


X] e~ st CT a ^(S, 0) = 1. 

o:—a_|_ ,ot — 


( 114 ) 














